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ABSTRACT 

This  report  is  the  first  in  a series  of  five  describing  an 
extensive  series  of  model  tests  on  air  blast  propagation 
in  underground  ammunition  storage  sites.  A -general  descrip- 
tion is  given  of  the  scope  and  purpose  of  the  tests  as  well 
as  a presentation  of  the  principles  of  scaling  and  various 
intrinsic  uncertainties  which  may  cause  problems  in  the 
modelling  technique.  The  report  also  surveys  the  test 
programme,  instrumentation,  and  data  reduction  for  the  whole 
test  series. 
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I . INTRODUCTION 

This  report  is  the  first  in  a series  of  five  / 1-6/* 
describing  an  extenive  series  of  model  tests  on  air 
blast  propagation  in  underground  ammunition  storage 
sites , 

These  tests  are  a natural  continuation  of  earlier 
work  performed  by  this  office  in  this  general  field 
/7-12/. 

I I Study  Objectives 

The  basic  objectives  for  performing  these  tests  have 
been  to  use  several  empirical  and  theoretical  scaling 
techniques  to  combine' air  blast  data  into  a form 
that  allows  the  prediction  of  blast  wave  propagation 
inside  a wide  range  of  typical  underground  ammunition 
storage  sites  as  a function  of  explosive  quantities 
and  geometrical  parameters.  These  methods  of  pre- 
dictions will  thus  give  the  design  engineer  a tool 
to  evaluate  the  effects  of  different  layouts  and 
geometries . In  particular,  it  will  be  possible  to 
specify  the  necessary  strength  of  blast  doors,  blast 
valves  and  blast  walls  to  prevent  the  propagation  of 
detonation  from  one  chamber  to  another  in  a connected 
chamber  storage  site.  Furthermore,  in  the  determina- 
tion of  the  pressure  and  impulse  at  the  exit  of  the 
tunnel  leading  into  the  underground  storage  complex, 
it  will  be  possible  to  determine  the  hazardous  area 
around  such  a site  /13,  14/. 

1 . 2 Problem  Background 

The  detailed  analysis  of  accidental  explosions  in 
underground  ammunition  storage  sites  in  regard  to 
blast  v/ave  propagation  generally  involves  three 
distinct  steps: 

Determination  of  the  pressure-time  history  in  the 
storage  chamber;  determination  of  the  blast  wave 
propagation  from  the  storage  chamber  down  the  main 
passageway  (or  between  chambers  in  connected  chamber 

* Numerals  in  / / designate  appended  references 
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storage  sites);  and  finally,  mapping  of  the  blast 
wave  emerging  from  the  exit  of  the  tunnel  in  the 
terrain  outside.  'This  analysis  can  generally  be 
performed  in  three  ways:  full  scale  testing;  theoreti- 

cal calculations;  or  finally,  model  testing. 

As  far  as  theoretical  treatment  is  concerned,  one 
can  in  principle  describe  the  blast  propagation  by 
sets  of  mathematical  equations,  and  attempt  to  solve 
these  equations  numerically.  However,  the  methods  of 
computation  are  often  so  involved  that  one  cannot 
economically  repeat  these  computations  while  varying  in 
a systematic  manner  all  of  the  physical  parameters 
which  may  affect  the  blast  wave.  One  important  point 
to  remember  is  that  one  has  no  idea  of  the  accuracy 
of  such  theoretical  solutions  until  proof  tested. 

The  second  possibility  is  experimental  studies  in 
large  or  full  scale.  Such  tests  would,  however,  be 
very  costly,  take  long  time  and  also  cause  instru- 
mentation and  manning  problems.  This  would  limit 
the  number  of  tests  to  a very  f ew  special  cases,  and 
there  could  be  a tendency  to  draw  conclusions  outside 
the  range  of  the  actual  experimental  and  geometrical 
conditions. 

The  third  and  last  possibility  is  model  testing 
using  similarity  methods.  This  method  forms  the  basis 
for  the  experiments  presented  in  the  present  reports. 

Scope  of  Report 

The  general  principles  of  scaling  are  presented  in 
Sec. 2 together  with  a discussion  of  various  intrinsic 
uncertainties  which  may  cause  problems  in  the  modelling 
technique.  The  test  plan  and  models  are  described 
in  Sec. 3 followed  by  a presentation  of  the  instrumenta- 
tion in  Sec. 4.  Section  5 contains  a description  of 
the  data  reduction  and  a discussion  of  the  experi- 
mental errors . 
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The  appendices  provide  more  detailed  sources  of  in- 
formation of  the  data  reduction  and  analysis  for 
this  whole  series  of  reports. 

Because  several  terms  that  relate  to  Blast  propaga- 
tion in  underground  ammunition  storage  are  used 
differently  by  different  authors,  a list  of  defini- 
tions is  also  included  as  a separate  appendix. 

2.  ■ THEORETICAL  CONSIDERATIONS 

2 • 1 Scaling  Relationships 

The  basis  for  the  scaling  relationships  used  in  the  ■ 
present  work  correspond  to  the  familiar  Hoplcinson 
scaling  laws  / 1 5 /.  The  use  of  models  requires  com- 
plete geometrical  similarity  between  the  model  and  the 
full  scale  system.  If  the  linear  scale  factor  is  n, 
the  relationship  between  distances  in  the  full  scale 
(subeript  E)  and  model  (subscript  M)  is 


LP  “ nIjM 

Likewise,  areas  scale  as 


(2.1  a) 


A — n A , 


(2.1  b) 


and  volumes  as 


Vp  = n 


(2.1  c) 


Eor  the  energy  release 


Ej,  = n^EM 


(2.1  d) 


This  energy  scaling  is  also  valid  for  the  explo- 
sive quantity  if  the  energy  release  per  unit  weight 
is  constant: 


Q1?  - n Qm 


(2.1  e) 


The  "blast  wave  parameters  scale  as 

Pp  ^ PM  (2.1  f) 

for  peak  pressure,  and 

I p = ( 2 . I g ) 

for  impulse,  and 

"kp  = (2.1  h ) 

for  characteristic  times. 

Inherent  in  this  simple  scaling  is  the  assumption 
of  perfect  gas  behaviour  of  air  and  no  differences 
in  ambient  conditions  (e.g.  temperature  and  pressure) 

' between  the  model  and  prototype.  Furthermore,  the 
theory  does  not  account  for  energy  dissipating  effects 
due  to  tunnel  or  chamber  walls  such  as:  thermal  energy 

transmission  to  the  walls;  elastic  or  inelastic  de- 
formation of  the  walls;  viscous  loss  due  to  wall  rough 
ness . 

If  one  tried  to  scale  all  these  effects  properly  from 
a full  scale  system,  it  would  require  tests  in  scale 
1:1.  It  is  therefore  important  to  try  to  put  an  upper 
limit  on  the  uncertainties  involved  in  neglecting 
non-scaled  effects.  This  will  be  attempted  in  Sec. 2.2 
2.4,  and  it  will  be  shown  that  the  effects  from  these 
sources  will  he  relatively  unimportant  except  maybe 
for  the  viscous  loss  due  to  wall  roughness.  However, 
as  the  relative  wall  roughness  is  normally  significantl 
larger  in  a full. scale  site  than  in  a steel  model,  the 
model  data  would  tend  to  he  conservative  or  on  the  safe 
side . 

As  will  be  discussed  in  the  succeeding  reports,  care- 
ful inspection  of  the  experimental  results  from  the 
various  model  tests  have  produced  extended  empirical 
scaling  relationships  representing  a wide  range  of 
important  geometrical  parameters.  One  should  note 
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the  distinction  between  scaling  laws  and  scaling 
relationships  in  this  connection  as  the  former  re- 
quire complete  geometrical  similarity  between  the 
model  and  the  prototype. 

2 • 2 Effects  of  Thermal  Interactions  with  the  Walls 

The  detonation  of  explosives  in  confined  regions  such 
as  storage  chambers  produces  hot  gases  from  which 
energy  is  transferred  by  heat  flow  to  the  much  cooler 
confining  walls.  There  are  two  mechanisms  involved 
in  this  heat  transfer:  Radiation  and  a combination 

of  -conduction  and  convection.  In  that  the  surrounding 
medium  is  rock  in  the  full  scale  case  and  steel  in  the 
present  model  tests, the  heat  losses  do  not  scale  due 
to  the  large  difference  in  thermal  properties. 

The  question  is  therefore  whether  these  effects  will 
be  of  such  magnitude  that  it  will  introduce  large 
systematic  errors  in  the  scaled  data. 

Earlier  theoretical  calculations  have  been  carried  out 
to  estimate  the  heat  conduction  in  a detonation 
chamber  in  full  scale  (rock)  and  steel  models  in 
linear  scales  1:10  and  1:100  / 1 6/ . The  basic  results 
are  reproduced  in  Pig. 2. 2 for  a loading  density  of 
100  kg  TNT/mb . It  can  be  seen  that  for  the  full 
scale  case  (chamber  radius  R = 6 m),  the  heat  loss 
is  relatively  unimportant.  Por  a steel  model  in  a 
linear  scale  1:100,  the  heat  loss  is  estimated  to  be 
about  14$  of  the  total  heat  at  a scaled  time  t/R.  = 

V 

20  ms/m,  corresponding  to  120  ms  in  the  full  scale 
case  / 1 7/ » Qualitative  arguments  have  also  been 
given  to  the  effect  that  the  relative  heat  loss  in- 
creases with  decreasing  loading  densities  and  that 
the  result  from  e.  1:100  steel  model  at  low  loading 
densities  are  expected  to  be  quite  misleading  unless 
a heat  loss  correction  is  applied,  fortunately,  the 
present  series  of  experiments  offers  some  results 
which  can  be  used  to  check  these  conclusions.  How- 
ever, first  an  attempt  will  be  made  to  convert  the 
thermal  loss  into  pressure  effects.  This  cannot  be 
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performed  in  a straightforward  way,  but  to  first 
order  the  pressure  decrease  in  an  invented  chamber 
due  to  conduction  loss  can  be  roughly  described  as 
/IS/ 


dp 

dt 


.i)/v. 


(2.2  a) 


Here,  dE  /dt  is  the  heat  loss  in  watts,  V is 
the  chamber  volume,  and  y = c /c^  the  usual 
specific  heat  ratio  for  the  detonation  products. 
Assuming  ^ to  be  constant,  the  relative  pressure 
decrease  can  be  found  from  Eq.  (2.2  a)  to  be 
approximately  equal  to  the  relative  heat  loss 


<*£  rO 
P 


(2.2  b) 


/mother  heat  transfer  which  produces  pressure  . 
decay,  is  radiation.  There  are  indications  that 
conduction  effects  at  the  high  temperatures  of 
internal  explosions  are  less  than  those  of  radia- 
tion /18/.  However,  the  radiation  loss  is  expected 
to  produce  approximately  the  same  effects  whether 
the  confining  walls  are  rock  or  steel. 

Relatively  good  experimental  estimates  of  the 
maximum  effects  on  peal:  pressure  and  positive 
impulse  of  the  heat  loss  are  provided  by  the 
closed  bomb  tests  in  Report  II  A / 1 / , and  the 
connected  chamber  storage  tests  in  Report 
IV  A /2/.  Details  are  given  in  these  reports 
and  only  one  example  will  be  discussed  here.  In  one 
typical  closed  bomb  test  with  a loading  density  of 
11  kg  INf/m3  (effective  chamber  radius  R Q — 

0,15  m),  it  was  found  that  the  pressure  fell 
approximately  15 7°  in  the  first  t/R  = 20  ms/m 
scaled  time,.  In  the  following  t/R  = 200  ms/m 
the  pressure  fell  only  about  5m. 
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These  pressure  drops  are  considered  to  be  upper 
limits  for  the  effects  from  heat  loss  due  to  con- 
duction in  as  much  as  part  of  the  pressure  drops 
found  in  these  tests  are  in  fact  due  to  small  leaks 
in  the  detonation  chamber  (Sec.3.5)as  well  as  radia- 
tion loss . As  noted  earlier,  the  effects  from 
radiation  probably  exceed  those  of  conduction.  Prom 
this  particular  test  and  from  similar  observations 
of  the  other  test  results  referred  to  earlier,  it  is 
therefore  concluded  that  the  heat  loss  will  be  of  no 
practical  importance  as  to  the  correctness  of  scaling 
the  present  steel  model  results  to  full  scale  in- 
stallations in  rock.  Thus,  the  calculations  in 
Ref.  IS  appear  to  overestimate  the  effects  from  heat 
loss  in  this  connection. 

One  possible  explanation  is  that  the  steel  walls  are 
assumed  to  he  completely  clean.  This  is  far  from 
being  true  under  normal  experimental  conditions.  The 
walls  will  be  somewhat  corroded  and  covered  with 
soot  due  to  the  incomplete  detonation  of  TUT.  This 
will  certainly  reduce  the  effective  heat  transfer  to 
the  steel,  hut  it  is  hardly  possible  to  be  more 
spesific  on  this  point  without  detailed  calculations. 

2 * 5 Mechanical  Interactions  with  the  Walls 

It  is  of  interest  to  put  some  upper  limits  on  the 
effects  from  elastic  or  plastic  deformations  of  the 
chamber  walls  in  the  present  tests.  Steel  is  elastic 
up  to  approximately  200  kbar  and  there  will  be 
negligible  direct  energy  loss  at  the  present  pressure 
levels  (Report  II  A).  However,  due  to  the  oscilla- 
tions of  the  pressure-time  history  (Report  II  A), 
the  chamber  -walls  will  also  oscillate.  This  will 
cause  some  changes  in  the  blast  wave  emerging  from 
the  chamber  with  a amah.-  lowering  of  the  pressure 
and  an  increase  in  the  duration  of  the  blast  wave , 
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For  the  simplified  case  of  a spherical  chamber  shown 
in  Fig.  2.3,  it  can  he  shown  that  the  deformation 
energy  from  a step  load  is  approximately  given  .by 


(2.3) 


with  parameters  defined  in  Fig.  2.3. 

To  evaluate  E^  relative  to  the  total  energy  release, 
one  particular  case  will  he  considered  with  the 
following  values  / 1 9/ : 


Q = 51,5  g TNT 
-< 

V ^ = 300  cm  = volume  of  detonation  chamber 
¥2  = 24500  cm  = total  volume  of  chamber 

This  produces 

Q/V^  — 172  kg/m3  - loading  density 
with  an  average  chamber  pressure  of  (Report  II  A) 
p — 2000  bar 

The  total  energy  release  is 

Ef.  rJ  1,4  x 105  Joule  (2700  Joule/g  TNT). 

The  volumes  and  correspond  to  the  following 
equivalent  spherical  radii: 

R1  = 4,2  cm 

R,^  = 18,0  cm . 

Using  these  values  and  the  values  given  in  Fig.  2.3 
for  steel,  Eq.  (2.3)  produces  the  total  deformation 
energy 


E 


D 


58  Joule 
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and  relative  to  the  total  energy  release 
• ED/Stre  4 x 10"4  or  0,04?$. 

As  p ocE^this  shows  that  the  mechanical  interaction 
with  the  walls  is  expected  to  have  a negligible 
effect  on  the  front  pressure. 

This  can  be  compared  with  a corresponding  full  scale 
chamber  in  granite  (with  the  same  geometry  as  in 
Pig.  2.3)  in  a linear  scale  100:1  of  the  model.  Using 
the  material  constants  in  Pig.  2.3,  the  relative  de- 
formation energy  is  approximately 
E j j/  E_j_  r\i  0,4  ia 

or  ten  times  as  high  as  for  the  steel  model.  (This 
can  be  seen  immediately  as  the  elastisity  module 
for  steel  is  approximately  ten  times  as  high  as  the 
elastisity  module  for  granite). 

The  general  conclusion  which  can  be  drawn  from  these 
approximate  calculations  is  that  the  energy  loss 
due  to  elastic  deformation  of  the  chamber  walls  will 
not  significantly  affect  the  scaling  relationships 
between  the  present  model  tests  and  prototype  in- 
stallations in  rock. 

The  proceeding  analysis  is  based  on  the  assumption 
of  elastic  deformations  of  the  detonation  chamber. 

In  practice,  however,  there  will  also  be  small 
permanent  deformations  of  the  chamber  and  exit  tube 
(visual  inspection  after  many  shots).  It  is  reason- 
able to  assume  that  this  will  produce  negligible 
energy  losses.  This  will  also  be  the  case  for  proto- 
type installations  in  rock  with  sufficient  overhead 
cover  to  prevent  venting,  as  the  destruction  is 
expected  to  take  place  in  the  vi cinity  of  cracks  and 
flaws  in  tiie  rock.  Even  with  the  destruction  of 
inner  concrete  liner  walls,  the  relative  energy  loss 
is.  expected  to  be  less  than  1?$  / 20/. 


In  the  model  tests,  shock  waves  in  the  air 
surrounding  the  model  will  also  he  generated.  The 
energy  of  .these  must  he  transmitted  through  the  steel 
walls,  and  has  to  he  less  than  the  deformation  energy 
taken  up  hy  the  model.  This  type  of  energy  loss  can 
therefore  also  he  neglected. 

2'. 4 Viscous  Losses  to  the  Walls 

There  are  numerous  experimental  results  pointing  to 
the  effect  that  blast  wave s in  tubes  and  tunnels 
may  he  damped  significantly  due  to  viscous  losses 
stemming  from  wall  roughness  / 21,  22/.  In  particular, 
the  results  from  recent  large  scale  tests  in  rock 
tunnels  in  Sweden  / 23/  and  Norway  /24/  have  showed  a 
viscous  damping  of  typically  40-80$  in  the  average 
peak  overpressure  for  travel  distances. of  about  20 
tunnel  diameters.  These  unusually  large  damping 
effects  were  due  to  large  wall  roughnesses  which 
were  typically  5-10$  of  the  tunnel  diameter. 

, For  the  steel  models  in  the  present  experiments,  the 
viscous  damping  will  he  much  less  as  the  tube  walls 
are  relatively  smooth  (Sec.  5).  This  has  been 
clearly  demonstrated  in  earlier  tests  of  one- dimens iona 
blast  wave  propagation  in  tubes  and  tunnels  with  dia- 
meters ranging  from  5 to  600  cm  / 21/.  Considering 
the  practical  accuracies  aimed  at  in  these  types  of 
experiments,  it  was  found  that  the  viscous  damping 
could  be  described  adequately  using  a model  originally 
proposed  by  Porzel  / 25,  22/.  This  model  has  also 
been  relatively  successful  in  describing  the  viscous 
damping  in  the  Swedish  tests  referred  to  earlier  /2o /. 
It  is  therefore  of  interest  to  outline  the  results 
of  this  model. 

From  thermodynamics  it  can  be  shown  that  pressure  re- 
duction due  to  wall  roughness  is  given  to  a good 
approximation  by  /2 5/: 
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Y(p)  = constant  - ^ 

Here,  Y(p)  is  an  impedance  function 

Y(p)  = lnp  -^7'T  ln  (p  +/"+  r)  ' -f 

where  p is  the  dimensionless  overpressure  ratio 

P = P^/p^  and  p^  = ambient  pressure  /27/.  Further- 
so  o o 

more,yois  a thermodynamic  parameter 

/*«  1) 

where  ^ is  the  usual  ratio  of  the  specific,  heats, 

V = cp^cv'  cons‘*:an't£  in  Eq.  (2.4  a)  is  a 

proportionality  factor  or  efficiency  with 

1/2  ^ 6 - 1 


It  has  been  argued  that 6 £2  1 for  "strong  impedance", 
i.e.  in  early  shock  wave  formation  and  6 ei  1/2  for 
"weak  impedance"  or  beyond  a so-called  choke  forma- 
tion in  the  tunnel.  e in.  Eq.(2.4  a)  is  the  average 
roughness  of  the  tunnel.  For  engineering  purposes, 
an  adequate  value  can  be  found  from 


e 


(Area)/ 


(Area) 


v/here  e equals  the  depth  of  roughness  and  the  area 
is  integrated  over  the  wall  of  the  tunnel.  Finally, 
D = hydraulic  diameter  = 4 A/S,  with  A the  cross- 
sectional  area  and  S the  tunnel  perimeter. 

For  high  pressures,  p >>  yA,(  = 6 for  air  with  y = 1,4 
Eq.(2.4  a)  reduces  to  Y(p)  ---  lnp  and  therefore  inp  = 

6g  L 

* p , which  has  the  solution 


(2.4  a) 
(2.4  b) 


(2.4  c) 


(2.4  d) 


(2.4  e) 


constant 
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L - L 

p(L)  = p(L0)exp(-k  — p — °)  (2.4  i) 

with  an  attenuation  constant  /28,  29/ 

^ = " (2.4  g) 

It  is  also  interesting  to  note  for  low  pressures, 
p«  Y ^ - i/p,  the  relative  pressure  attenuation 
expressed  in  Bq.(2„4  a)  is  minutely  small  compared  to 
the  high  pressure  attenuation  in  Eq.(2.4  f). 

The  value  for  e in  Eq.(2.4  a)  has  been  fitted  to  experi- 
mental data  referred  to  earlier  / 21,  30/  and  the  results 
show  fair  agreement  with  the  values  for  e when  measured 
directly. 

2 • 5 Validity  of  Model  and  Large  Scale  Testing 

It  is  often  argued  that  one  has  no  check  of  the 
accuracy  of  model  tests  unless  a direct  model- 
prototype  comparison  is  conducted „ In  as  much  as  there 
are  relatively  few  large  scale  tests  on  underground 
ammunition  storage,  only  a few  comparisons  for  very 
special  cases  can  be  made.  It  will  therefore  be  con- 
sidered very  important,  throughout  this  series  of 
reports  to  discuss  in  some  detail  those  comparisons 
for  which  both  model  and  prototype  data  are  available. 

The  term  "validity  of  model  testing"  is  often  misused 
as  one  receives  the  impression  that  model-prototype 
comparison  is  always  needed  because  the  analysis 
applies  similtude  technique.  One  should  l’ather  con- 
sider the  few  model-prototype  comparisons  as  tests 
of  the  assumptions,  and  accept  the  validity  of  model 
tests  with  similar  accuracy  for  analogous  problems. 


Another  important  point  to  remember  is  that  the  model 
tests  described  in  these  reports  apply  methods  of 
obtaining  engineering  answers  to  questions  which  cannot 
be  found  more  accurately  at  an  acceptable  cost  any 
other  way.  The  very  few  large  scale  tests  conducted 
up  to  the  present  time  are  not  sufficient  by 
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themselves  to  form  a basis  for  safety  criteria,,  There 
are  several  reasons  for  this: 

- There  is  a lack  of  systematic  investigations  for  differ- 
ent geometries 

- The  majority  of  the  test  sites  are  not  geometrically 
similar  to  sites  in  use 

- Energy  losses  due  to  friction  and  viscous  interaction 
with  the  walls  are  expected  to  be  significantly  different 
for  several  of  the  tests  due  to  relatively  large  differ- 
ences in  wall  roughness  from  site  to  site.  This  compli- 
cates the  quantitative  comparisons  considerably 

- There  are  severe  difficulties  connected  with  the  instru- 
mentation due  to  the  destructive  power  from  debris  and 
the  need  for  long  transmission  lines  between  the  pressure 
transducers  and  the  recording  instruments.  The  inter- 
pretation of  many  of  the  test  results  has  therefore  been, 
quite  uncertain 

- Attempts  to  confirm  the  test  results  theoretically 
have  so  far  been  inconclusive  as  to  the  possibility  of 
extending  the  results  to  widely  different  geometries  and 
providing  reliable  information  of  the  uncertainties  in- 
volved 

From  these  general  observations  it  therefore  appears  that 
the  sparse  data  from  the  large  scale  results  should  be 
used  with  some  caution.  This  is  especially  true  in 
attempting  to  generalize  or  extend  the  results  to  a wide 
range  of  different  geometries.  These  tests  will,  how- 
ever, be  quite  important  fox-  direct  comparisons  with 
geometrically  similar  model  tests.  In  the  event  of  future 
costly  large  scale  experiments  it  is  considered  essential 
to  have  these  preceded  by  model  tests  as  an  aid  in  the 
planning. 
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3. 


■TEST  PROGRAMME  AND  MODELS 

The  test  programme  was  divided  into  four  major  parts 
which  are  reported  in  the  four  remaining  reports  of  this 
series: 


3*1  ?©P2£j..II 

Closed  bomb  tests  to  determine  the  average  pressure  exerted 
on  the  chamber  walls  in  underground  ammunition  storage 
sites.  The  parameters  which  were  varied  in  these  tests 
are  shown  in  Fig.  3. It 

V-  = chamber  volume 

Q.  = charge  weight  and  type  of  explosive 

A- = vented  area 
J — 

3.2 

De termination  of  the  blast  wave  propagation  in  the  passage- 
way of  a single  chamber  storage  site.  The  parameters 
which  were  varied  in  these  tests  are  shown  in  Fig.  3.2: 

V^=  chamber  volume 

Q = explosive  charge  weight  (TNT) 

A.=  passageway  cross  section 
J 

L = distance  along  passageway 


3.3  Report_IV 

Determination  of  .the  blast  wave  propagation  in  the  main 
passageway  of  a connected  chamber  storage  site.  The  para- 
meters which  were  varied  in  these  tests  are  shown  in 
Fig.  3*3: 

= chamber  volume 

Q - explosive  charge  weight  (TNT) 

I 

A.  = branch  passageway  cross  section 

J 

A,  = main  passageway  cross  section 
vm  - angle  between  branch  and  main  pas  sage 'way 
- distance  along  main  passageway 


1 
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3.4  Report^V 

Determination  of  the  blast  wave  propagation  between  connected 
chamber  storage  sites  and  the  blast  load  on  doors  in  three 
sites.  The  different  geometries  and  the  parameters  which 
were  varied  in  the-  tests  are  shorn  in  Fig.  3.4. 

Q = explosive  charge  weight 
I = distance  along  main  passageway 

3 . 5 Model  Construction 

The  chambers  were  the  most  critical  parts  in  the  construc- 
tion of  the  models  in  regard  to  strength  as  average 
pressures  up  to  200  'bars  were  expected  with  pressure  peaks 
up  to  8 times  this  value.  This  required  a type  of  steel 
which  would  not  be  too  brittle  and  a quality  St  50  v/as 
chosen.  Example  of  the  design  of  one  particular  chamber 
. (V  = 15200  cm"j  is  shown  in  Fig.  3.5  a with  the  standard 
codes  used  in  Norway. 

To  vary  the  chamber  volume,,  various  insertions  were  used  and 
some  typical  examples  are  shown  in  Fig.  3.5  b.  As  may 
be  seen,  the  inlays  were  mounted  directly  on  to  one  of 
the  endpieces  of  the  chamber. 

Examples  of  the  -connections  between’  chambers  and  the  tubes 
are  shown  in  Fig.  3.5  c.  This  figure  also  shows  the  design 
of  some  of  the  straight  tubes  (main  passageway)  with  the 
pressure  gauge  mountings. 

All  chambers,  branch  passageways  and  main  passageways 
were  cylindrical,  but  using  the  cross  section  of  the  main 
passageway  as  the  scale  reference,  these  steel  models  are 
ideally  thought  to  be  in  a linear  scale  from  1:50  to 
1:100  of  typical  configuration  of  Norwegian  underground 
storage  sites.  The  difference  in  the  geometry  of  the 
cross  sections  between  model  and  full  scale,  is  thought 
to  have  minor  influence  upon  the  test  result. 
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4.  INSTRUMENTATION 

4 . 1 Pressure  'Transducers  and  Measuring  Chain 

To  measure  the  pressure- time  history  in  the'’ various  models, 
standard  instrumentation  techniques  were  employed.  The 
pressure  transducers  in  the  detonation  chambers  and  exit 
tubes  were  of  the  Iiistler  603  B,  601  H,.  and  6201  type  with  rise 
times  of  100  microseconds  and  a natural  frequency  of  400  kHz. 

The  mounting  of  the  transducer  in  the  detonation  chamber, 
as  shown  in  Pig.  4.1  a,  was  especially  designed  to  reduce 
the  thermal  response  of  the  transducer  to  the  high 
temperatures  of  the  detonation  products.  To  achieve  this, 
a relatively  large  cavity  in  front  of  the  transducer  was 
filled  with  a mixture  of  silicone  oil  and  graphite.  This 
arrangement  is  expected  to  reduce  the  rise  time  of  the 
transducer  by  a-  factor  of  about  30,  but  this  will  not  be 
critical  for  the  type  of  experiment  reported  here. 

Similar  precautions  were  taken  in  mounting  the  transducers 
in  the  tubes  as  shown  in  Pig.  4.1  b,  where  a 2 mm  thick 
Silicon  rubber  plate  served  as  a heat  shield.  This  pro- 
tection proved  to  be  sufficient  as  the  temperature  rise 
v/as  significantly  smaller  in  the  tubes  than  in  the  chamber. 

The  signals  from  the  pressure  transducers  were  processed 
as  shown  in  Fig.  4.1  c.  The  output  from  the  charge  ampli- 
fiers were  recorded  in  either  of  two  different  modes: 

~ Mode  A:  The  pressure- time  histories  were  photographed 
directly  from  oscilloscope  displays 

- Mode  B:  The  pressure  recordings  were  stored  on  an  analog 
tape  recorder  and  subsequently  played  back  to  a CSC  U.V, 
writer  using  light  sensitive  paper. 

Mode  A gives  intrinsically  the  best  resolution  of  about 
5 sec,  compared  to  an  approximate  resolution  of  25^  sec 
for  mode  B (40  kHz  bandwidth  of  Ampex). 


4.2  Transducer  Calibration  and  System  Linearitv 

— - i — - - _ . . ------ - — ■■■-■'-  ■■  ■■ 

The  calibration  data  for  the  transducers  supplied  by  the 
manufacturer  were  used  to  obtain  the  pressure.  Several 
transducers  were  occasionally  mounted  together  in  an  explo 
sive  charge  driven  shock  tube  to  check  the  overall  relativ 
accuracy  of  these  calibrations.  This  is  discussed  in 
Appendix  B. 

In  spite  of  relatively  reliable  calibration  data  on  the 
transducers,  there  will  of  course  be  errors  Introduced  by 
the  other  components  in  the  measuring  chains.  Some  idea 
of  these  uncertainties  can  be  estimated  from  the  specifica 
tions  given  in  Table. 4. 2. 

5.  DATA  REDUCTION 

5 • 1 Curve  Tracing  and  A/D  Conversion 

In  order  to  effectively  analyse  the  pressure-time  record- 
ings, a D-mac  curve- tracer  was  used  to  convert  the  analog 
signals  to  digital  form.  This  is  discussed  in  Appendix  C. 
Further  data,  reduction  was  performed  on  an  IBM  570/155  and 
frontpressure  (maximum  average  overpressure),  impulse, 
positive  duration  and  time  of  arrival  were  obtained.  Thes 
data  were  then  sorted  according  to  test  configuration, 
charge  weight  etc,  and  printed  in  tabular  form  / 6/  and 
punched  on  cards  for  further  analysis. 


5 • 2 Fitting  Technique  for  Scaling  Relationships 

For  the  27  different  model  configurations  described  in 
Reports  III  and  IV,  a.  total  of  2200  blast  wave  parameter 
values  (peak  pressure,  impulse,  and  positive  duration) 
were  obtained.  Attempts  were  made  to  fit  these  data  'jo 
several  equations  based  on  various  seeding  relationships 
to  link  the  data  from  these  different  co:.if i .mirations  to- 


gether. 
normal 


■i.  ..iVj  V' 


procedure 


13  type  oi  equations 

r 


tested  am  not  per/rre  tn c 
juo  of  least  squares  fitting  and  so  an  iierriit 
•/Rich  requires  no  derivatives  was  used.  This  is 


discussed  in  more  detail  in  Anoendix  D. 
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5 • 3 Electrical  Filter  for  Chamber  Pressure 

As  will  be  discussed  in  Report  II,.  the  recorded  signals  from 
the  transducers  mounted  in  the  detonation  chamber  contained 
high  frequency  ringing.  This  was  mainly  due  to  the  multiple 
reflected  wave s from  the  walls,  but  also  to  some  extent 
due  to  vibrations  of  the  chamber  walls  /31/.  The  signals 
from  the  magnetic  tape  were  in  this  case  passed  through 
a carefully  designed  electrical  filter  used  elsewhere 
/ 32/ . There  were,  however,  some,  alterations  to  the  filter 
used  in  Ref. 32  and  a full  description  is  therefore  given 
in  Appendix  E. 
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APPENDIX  A 


DEFINITIONS  OF  SYMBOLS 


Cross-sectional  area  defined  as: 

A.  main  passageway,  single  chamber  storage  site 

u 

I 

A..  branch  passageway,  connected  chamber  storage  si 

v 

A^  main  passageway,  connected  chamber  storage  site 


c Velocity  of  sound 

D = 4 A/S,  hydraulic  diameter. 

E Elasticity  module 

E Conduction  heat  loss 

c 

E-p  Deformation  energy 

E^  Total  energy  release  of  explosive 

e local  height  of  wall  roughness  element 
'<j  Average  height  of  wall  roughness  element 

I Impulse  (bar-  ms) 

k = 2 G e/D,  attenuation  constant 

1 Distance  along  main  passageway  from: 

a)  Chamber  exit  in  single  chamber  storage  sites 

b)  Exit  of  branch  passageway  in  connected  chamber 
storage  sites 


n Linear  scaling  factor 
p Peak  overpressure  (bar) 


Q TNT  charge  weight 

S Perimeter  of  cross-section  of  main  passageway 
t Positive  duration  of  pressure  pulse 

t Arrival  time  of  shock  front  from  time  of  detonation 

St 


A 2 


vm  Angle  ‘between  branch  and  main  passageway 

V Chamber  volume  (closed  bomb) 

V.  Chamber  volume  for  single  and  connected,  chamber  storage 
sites 

V,  Total  effective  volume  (Volume  of  chamber  and  tunnel 
system  up  to  the  observation  point 

X Energy  division  factor  for  connected  chamber  storage 

sites 

Y Impedance  function 


= C /Cv  , specific  heat  ratio 
Efficiency  of  feedback  energy 
Lamp's  constant,  def.  in  Pig.  2.3 
Lamp's  constant,  def.  in  Eig.  2.3 

- ( v/+  o/(^-  o 


DEFINITION  OF  TERMS  USED  IN  UNDERGROUND  STORAGE 


Adit:  A passage  or  tunnel  leading  into  an  underground 

storage  complex. 

See  Underground  Storage 

After-Burning : Blast  or  flames  from  post-detonation  reaction 

HSetv/een  combustible  explosion  gases  and  ambient  oxygen. 
This  effect  occurs  for  oxygen  deficient  explosives  as 
for  example  T N T. 

See  Blast  Nave,  Jet-Effect 

Average  Overpressure : The  pressure/ time  curve  obtained  by 

Taking' ’Hie  mean  values  of  the  pressures  of  the  peaks 
and  troughs  in  a pressure  oscillating  system. 

Average  Peak  overpressure : The  maximum  value  of  the  average' 

o v e r p r ess u f e . Tee  "Average  Overpressure 


Ball  of  Fire  (or  Fireball):  The  luminous  sphere  o: 

G aTcTTT o rm s a few  millionths  of  a second  after 


hot  gases 
an  circle 


and  immediately  starts  to  expand  and  cool.  The  exterior 
the  ball  of  fire  is  initially  sharply  defined  by  the 
luminous  shock  front  (in  air)  and  later  by  the  limits  01 
the  hot  gases  themselves.  See  Breakaway ■ 


d/or 


Blast  Attenuation  : The  decrease  in  peak  overpressure  sir. 

impulse  versus  distance  from  the  explosion  source  due  • 
rarefaction  wave  and/or  due  to  energy  dissipating  ef fee 
in  tunnel  systems  such  as  viscous  loss  due  to  wall 


rouglmes; 


iherrnal  energy  transmission  to  the  walls,  ve: 


formation  etc. 

See  Rarefaction  Wave,  Viscous  Loss, 


.1  Roualuiess.  Ver{ 


Blast  Loading  : The  loading  (or  force)  on  an  object  caused  cp 
the  -air  "blast  from  an  explosion  striking  and  flowing 

a combination  of  overpressure 


around  the  object.  It  i 


(or  diffraction)  and  dynamic  pressure  (or  drag) 


See  Diffraction, 


Dr 


'yn  ami  c 


s sur 


Overrr 


Blast  protective  Construction 
sue : 


:-x>e  ci  f i call  v d e s i sn e cl  bu i 1 d i r 


element : 


ast  doors,  blast  valves 


c,  'co  re  as 

or  prevent  the  entry  into  a chamber  of  the  blast  wave  f 
the  explosion  of  the  contents  of  an  adjacent  chamber, 


See  Blast  Door 


: 1 V 


Blast 


Blast  Resist ins  Doi 


ioo: 


me  err 


specifically  des 
into  a chamber  of  the  *?».*? es*. 


signeo.  to  trover 

try  into  a cnamoer  oi  me  gases,  flame  and  air  sh: 
from  the  explosion  of  the  contents  of  an  adjacent  rdvr 


I 


Blast  Seal  inn:  laws  : Formulas  which  permit  the  calculation  of 

the  properties,  e.g. , overpressure,  dynamic  pressure,  time 
of  arrival,  duration,  etc, of  a blast  wave  at  any  distance 
from  an  explosion  of  specified  energy  from  the  "known 
variation  with  distance  of  these  properties  for  a reference 
explosion  of  known  energy,  e.g.  of  1 kilogram  TuT. 

See  Cube  Root  Law. 

Blast  Trap  i Special  construction  features  in  the  tunnel  system 
o f an  underground  storage  sits  such  as  turns,  crossovers, 
obstacles,  cross  section  changes  etc.,  to  reduce  the  blast 
wave  leaving  or  entering  a.  passageway. 

Blast  Valve  i A valve  in  the  ventilation,  system  specifically 
' designed  to  prevent  the  entry  into  a chamber  of  the  blast 
wave  from  the  explosion  of  the  content  of  an  adjacent  chamc- 

Blast  W all  : A wall  which  supports  (or  bears)  part  or  mass  of 

a structure,  such  as  a blast  door,  specifically  designed  to 
withstand  the  blast  wave  from  the  explosion  of  the  contents 
of  an  adjacent  chamber. 

Blast  Wave  : A pressure  pulse  of  air,  accompanied  by  winds, 
propagated  continuously  from  an  explosion. 

See  Shock  wave. 

Branch  Passageway  : The  tunnel  leading  from  the  storage  chamber 
to  tnc  main  passageway  in  a connected  chamber  storage  site. 

Breakaway  : The  onset  of  a condition  in  which  the  shock  front 

(in  the  air) moves  away  from  the  exterior  of  the  expanding 
ball  of  fire  produced  by  the  explosion. 

See  Ball  of  fire,  Shock  f rent . 

Bulk  Velocity  : The  resultant  velocity  imparted  to  a mass  by  a 
“ ."“series  'of  shock  and  reflected  waves. 

Burst  : Explosion  or  detonation. 


Chamber  Interval  : The  interval  between  the  natural  walxs 
, 'adTaceru "underground  storage  sites. 

See  Qu  an  t i t y fist  an  c s 

Chamber  Btornse  Site  : Normally  a chamber  excavated  in  c ' 
“ or  a. 'serf-.; s of  such  chambers  specially  adapted  for  v: 
storage  of  ammunition  and  explosives.  A cavern  may  e 
structurally  modified  to  be  used  as  a chamber  sv-craj 

Connected  Chamber  Storage  Site  \ An  excavated  c numbs r v/m 
^ ^ _ c 6nh  e cf  e ;i  "“oy"^aTi'“d  u c t s or  pasoagev;yys  to  anotn 

st o rage  chamber . 


Contained  Explosion  : An  explosion  in  an  underground  sto: 
grp EVTTEEl s or  no  venting  to  the  surface. 


Cover  : The  solid  ground  situated  between  the  ceiling  or  wall 
of  an  underground  chamber  and  the  nearest  exterior  surfac- 

Crater  i A hole  or  chasm  in  the  overburden  caused  by  a 
subterranean  explosion. 

Cube  Root  Law  : A scaling  law  applicable  to  many  blast 

■ phenomena.  It  relates  the  time  and  distance  at  which  a 
given  blast  effect  is  observed  to  the  cube  root  of  the 
energy  yield  of  the  explosion. 


Damage  Criteria  ! Standards  or  measures  used  in  estimating 
specific  levels  of  damage. 

Detonation,  A violent  chemical  reaction  within  a chemical 
compound  or  mechanical  mixture  evolving  heat  and  high 
pressures,  A detonation',  in  contradistinction  to  deflagra- 
tion, is  the  reaction  which  proceeds  through  the  reacted 
material  toward  the  unroacted  material  at  a supersonic 
velocity.  The  result  of  the  chemical  reaction  is  an  exer- 
tion of  extremely  high  pressures  on  the  surrounding  me-  ■ 
dium  forming  a propagating  shock  wave  which  is  originally 
of  supersonic  velocity.  A detonation,  when  the  material 
is  located  on  or  near  the  surface  of  the  ground,  is  nor- 
mally characterised  by  a crater. 

Diffraction  : The  bending  of  waves  around  the  edges  of  object: 
Th  connection  with  a blast  wave  impinging  on  a structure 
diffraction  refers  to  the  passage  around  and  envelopment 
of  the  structure  by  the  blast  wave.  Diffraction  loading  : 
the  force  (or  loading)  on  the  structure  during  the 
envelopment  process. 

Drag  Loading  : The  force  on  an  object  or  structure  due  to 
-'"the  "transient  winds  accompanying  the  passage  of  a blast 
The  drag  pressure  is  the  product  cf  the  dynamic  pressure 
and  a coefficient  which  is  dependent  upon  the  shape  (or 
geometry)  of  the  structure  or  object. 

See  Dyn  ami c o r e 3 s u r e . 


Dynamic  Pressure  : The  air  pressure  which  results  from  the 
mass' a rf~TTo w (or  wind)  behind  the  shock  front  of  a 
blast  wave.  It  is  equal  to  the  product  of  half  the 
density  of  the  a.ir  through  which,  the  blast  v;ave  passes 
and  the  square  of  the  particle  (or  wind)  velocity  in  the 
wave  as  it  impinges  on  the  object  or  structure. 


Effective  Net  Explosive  Quantity:  See  Net  Explosive  Quant  it 


Elastic  Range  : The  stress  range  in  which,  a material  will, 
recover  its  original  form  when  the  force. (or  loading) 
is  removed.  Elastic  deformation  refers  to  dimensional 
changes  occurring  within  the  elastic  range. 

See  Plastic  Range. 

Emergent  Shock  or  Blast  : The  shock  wave  issuing  from  an 
' cicTi t , a”  vent  or  a~'c rater. 

Explosion.  The  rapid  production  of  hot  gases  at  a high 
pressure  as  the  result  of  a chemical  reaction  and  the 
suddc-n  release  of  this  energy  to  cause  strong  dynamic 
stresses  in  the  surroundings.  The  term  usually  relates  to 
the  effect  of  a detonation  of  initiating  explosives  and 
high  explosives  out  applies  in  certain  circumstances  to 
,the  effect  of  a deflagration  of  -propellant  explosives. 

kjiP.b  Q s i on  Pro  pa  gat  i o n : A process  where  the  explosion  in 
one  chamber  causes  an  explosion  in  an  adjacent  chamber 
due  to  blast,  hot  flames,  or  debris. 

Explosive . A substance  manufactured  with  a view  to  produce 
a practical  effect  by  explosion  or  a pyrotechnic  effect. 

An  explosive  atmosphere  of  gas,  vapor  or  dust  is  not  con- 
sidered to  be  an  exp 3. o save . 


Ext  e r i o r Qu ant 1 1 y~ D i a t 


The  minimum  permissible 


distance  between  a building  or  stack  containing  ammunition 


or  explosives  and  a building  or  an  assembly 
public  traffic  route  outside,  the  explosive 


y pxace  or 


:a. 


Free  Air  Overpressure  (or  Free 
pressure ~Tn  exce  3 so" rTT 
created  in  the  air  by  the 


Air  Pressure)  : The  unre fleets 
Te~ am biewm  "Trnuo  spheri c pre s sure 
blast  wave  from  an  explosion. 


Fu  1 1 s c al  0 T r i a Is 


elves  (tons) 
simulate  and 


explosions . 


: Detonation  of  large  quantities  of  explo- 

in  large  scale  installations  in  rock  to 
observe  the  effects  from  possible  accidental 
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.Geometrical  Copy  Model  Tests:  Tests  in  a scaled  down  geometri- 

cally similar  rSoHel  of  a 'full  scale  installation.  The  linear 
scales  of  the  model  are  reduced  by  a scaling  factor  ns  and 
the  explosive  charge  weight  by  a factor  of  n3  compared  to 
the  full  scale  installation. 


Impulse : The  product  of  the  overpressure  (or  dynamic  pressure) 

from  the  blast  v/ave  of  an  explosion  and  the  time  during 
• which  it  acts  at  a given  point.  More  specifically,  iWis  rue 
integral,  with  respect  to  time,  of  the  overpressure  (or 
dynamic  pressure),  the  integration  being  be Ween  the  tine 
of  arrival  ox  the  blast  wave  and  that  at  which  the  over- 
pressure (or  dynamic  pressure)  returns  to  zero  at  the 
given  point. 


In t e pi  o r Q,u ant  1 1 ?.r- 1) i s t an c e . 
distance-  between  a building 


The  minimum  permissible 
or  stack  containing  ammuni- 


tion or  explosives  and  another  such  building  or  stack 
inside  the  explosives  area. 


]' , i n g sr  a c a 1 i n g Facto r : 
of  'aTuIT  scale  "Inst 
model.  The  model  te 
scale  1 ;n. 


The  ratio  n between  the  linear  dimension 
allation  and  a geometrically  similar 
sts  are  then  said  to  be  in  a linear 


Loading  : The  force  on  an  object  or  structure  or  ele 
a structure.  The  loading  due  to  blast  is  eaual 
net  pressure  in  excess  of  the  ambient  value  mul 
by  the  area  of  the  loaded  object  etc. 


mens 
to  t 
tipi 


Loading  density  : The  ratio  Q/V,  where  Q is  the  total  net 

explosive"  quantity  ana  V .is  the  volume  of  the  confine- 
ment, e.g,  a storage  chamber. 


Mach  front  : See  Mach  Stem. 


inch  Peri  or  : The  region  on  the  surface  at  which  1 me  ..non 
luisW’ormed  as  the  result  of  a particular  explosion 
air. 
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Kach  St  era  : The  shock  front  formed  by  the 

'rents 


tuoion  O-  xhs  arc 
om  an  explosion.  The  teem. 


and  reflected  shock 
generally  used  with  reference  to  a blast  wave,  propa- 
gated in  the  air,  reflected  at  the  surface  of  the  earth 
The  Mach  stem  is  nearly  perpendicular  to  the  reflecting 
surface  and  presents  a slightly  convex  (forward)  front. 
The  Kach  stern  is  also  called  the  Mach  front. 

See  Shock  front,  Shock  wave. 


Main  Pas sage wav  : A passage  or  tunnel  leading  into  an  under- 

“*  gr ouncPsfbrag e complex. 

Mass  Explosion.  An  explosion  which  affects  almost  in- 
stantaneously tne  entire  quantity  of  explosives  under 
consideration.  The  term  usually  relates  to  detonation 
but  also  applies  to  deflagration  when  the  practical  ef- 
fects are  similar,  for  example  the  mass  deflagration  of 
propellant  under  very  strong  confinement  so  as  to  pro- 
duce a cursting  effect  and  a serious  hazard  from  debris. 

Model  Tests  : See  Geometrical  Copy  Model  Tests 


Mot  Explosives  Quantity.  The  net  explosives  quantity  is 
the  total  explosives  contents  of  ammunition  unless  it 
has  been  determined  that  the  effective  quantity  is  signi- 
ficantly different  from  the  actual  quantity.  It  does  not 
include  such  substances  as  white  phosphorous,  war  gases 
or  smoke  and  incendiary  compositions  unless  these  sub- 
stances contribute  significantly  to  the  dominant  hazard 
of  the  hazard  division  concerned. 


Overhead  Cove”  : See  Cover 


Overpressure  : The  transient 
bounds  rer  sauarc  inch 


pressure,  usually  expressed  in  bars 
exceeding  the  ambient  pressure, 


manifest el  in  the  shock  (or  blast)  wave  from  an  explosion 
The  variation  of  the  overpressure  with  time  depends  on 
the  energy  yield  of  the  explosion,  the  distance  from  the 
point  of  burst,  and  the  medium  in  which  the  explosive  is 
detonated.  The  peak  overpressure  is  the  maximum  value  of 
the  overpressure  at  a given  location  and  is  generally 
experienced  at  the  instant  the  shock  (or  blast)  wave 
reaches  that  location. 

See  Shock  wave. 


Particle  Velocity  : The  local  velocity  imparted  by  the 
t ransmlssion  of  a shock  or  a reflected  wave . 

Peak  Overpressure  : See  Overpressure . 

Plastic  Range  : The  stress  range  in  which  a material  will  not 

fail  when  subjected  to  the  action  of  a force,  but  v;ill  r.o' 
recover  completely,  so  that  a permanent  deformation  resui" 
when  the  force  is  removed,,  Plastic  defer:-. '.op  ion  refers  to 
dimensional  changes  occurring  within  the  plastic  range. 
See  ; SI a s t i c range. 

Positive  Phase  : See  Shock  V/ave. 


Propagation  by  Hot  Gases  i See  Explosion  Propagation. 
Propagation  by  Rock  Spall  ! See  Explosion  Propagation. 


Q,u ant 1 v, y - P i s t an c o . The  minimum  permissible  distance  be- 
tween a building  or  stack  containing  a given  quantity 
of  ammunition  or  explosives  and  any  other  exposed 
structure,  assembly  place  or  traffic  route.  It  is  based 
on  an  acceptable  risk  to  life  and  property  from  the  ef- 
fects of  a mass  fire  or  an  explosion. 


Rarefaction  v/ave  ; An  underpressure  pulse  produced  as  tiie 

blast  wave  emerges  from  the  open  end  of  a duct  and  travels 
back  into  the  duct  increasing  the  normal  pressure  decay. 

Reflection  Pactor  : The  ratio  of  the  total  (reflected) 

pressure  to  the  incident  pressure  when  a shock  (or  blast’ 
wave  travelling  in  one  medium  strikes  another. 

Reflected  Pressure  : The  total  pressure  which,  results 

instantaneously  at  the  surface  when  a shook  (or  blast) 
wave  travelling  in  one  medium  strikes  another  medium, 
e.g. , at  the  instant  when  the  front  of  a blast  wave  in 
air  strikes  the  surface  of  an  object  or  structure. 


Re: ver.be re tions  : The  pressure  oscillations  produced  in  a. 

confined  space  by  the  transit  of  shocks  and  reflected 
shocks  in  the  gases. 
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Scaling  Factor  : See  Linear  Scaling  Factor. 


Scaling  Lav/:  A mathematical  relationship  which  permits  the 
effects  of  an  explosion  of  given  energy  yield  to  he 
determined  as  a function  of  distance  from  the  explo- 
sion, provided  the  corresponding  effect  is  kncv;n  as 
a function  of  distance  for  a reference  explosion, 
e.g.  1 kg  TNT.  See  Blast  scaling  law,  Cube  . oot  lav/. 


Secondarv  Peaks 


The  highest  pressure  recorded  in  the  trans- 


mission of  each  of  a series  of  reflected  shocks  following 
the  initial  shock  and  its  primary  reflected  shock. 

Shock  Front  (or  Pressure  Front)  : The  fairlv  sharp  boundary 

between  the  pressure  dis t urbane e created  by  an  explosion 
(in  air,  water  or  earth.)  and  the  ambient  atmosphere,  v/axer 
or  earth  respectively.  It  constitutes  the  front  of  the  shock 
• (or  blast)  wave. 


Shock  Wave  : A continuouslv  propagated  pressure  pulse  (or 


wave 


• ) 


water  or  eerti 


"in  The  surrounding  medium  which  mav  be  air. 
initiated  bv 

explosion.  A shock  wave  in  air  is  generally  referre 


because  it  is  similar  to  (and  is  accompanied  by) 


of  a shock  (or 


the  expansion  of  the  hot  gases  produced  in  an 

to  as 

blast  wave 

strong,  but  transient,  winds.  The  duration 
blast)  wave  is  distinguished  by  two  phases.  First  there 
the  positive  (or  compression ) phase  during  which  the 
pressure  rises  ver-  sharo.lv  to  a vslie  that  is  higher  than 
ambient  and  then  decreases  rapidly  to  the  ambient  pressure 
The  duration  of  the  positive  phase  increases  and  the  maxim 
(peak)  pressure  decreases  with  increasing  distance  from  an 
explosion  of  given  energy  yield.  In' the  second  Phase,  the 
negative  ( or  suction)  phase,  the  pressure  falls  below 
ambient  and.  then  returns  to  the  ambient  value.  The 


' *.  . i j.  C!  w . 


of  the  negative  phase  is  approximately  conacan.  Anro.ig^ou , 
the  blast  wave  history  and  may  be  several  times  the  dxvr.x i y 
of  the  positive  phase.  Deviations  from  the  ambient  pressur 
during  the  negative  phase  are  never  large  and  they  docrer-s 
with  increasing  distance  from  the  explosion. 

See  : 0 v e r o r e s sure. 


Sin/rlo  Chamber  Storrre  Site  : An  excavated  chamber  which  has 

the  exterior  and  is  not  connects.. 


its  own  entrance  ire- 
air  ducts  or  po  ssnme'.vavs  to  anv  other  storage 


nance.: 


Site.  Any  cavern,  chamber,  building,  cell  or  stack  which 
contains  ammunition  and/or  explosives. 


Specific  Impulse  : The  impulse  per  unit  area. 
Static  Pr e s s ir a : See  Ov er  or  c s s ir  e . 


St  ora  me  Chamber  : An.  excavated  chamber  in  rock  used  for  storage 

of  explosives  and  ammunition. 


\))  CD 
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Thermal^ Energy;  The  energy  emitted  from  the  ball  0f  fire 
thermal  radiation.  The  total  amount  of  thermal  encr 
received  per  unit  area  at  a specified  distance  from 
nuclear  (or  atomic)  explosion  is  generally  expressed 
in  terms  of  calories  jjer  square  centimeter. 

TNT  Equivalent:  A measure  of  the  energy  retested  in  the 

detonation  of  ammunition  or  explosives  of  a given 
quantity  expressed  in  terns  of  the  quantity  of  TNT 
which  would  release  the  same  amount  of  energy  when 
exploded.  The  TNT  equivalent  is  usually  stated  in  ideo- 
grams or  tons.  The  basis  of  the  TNT  equivalence  is 
that  the  explosion  of  1 ton  of  TNT  releases  109  calories 
of  energy. 


Unconta in ed  Explosion  : An  explosion  taking  place  in  open  ai] 
or- on. a free  surface. 

Underrt  aid  Explosion  : The  explosion  of  explosives  or 

i .unit ion  with  its  centre  beneath  the  surface  of  the 
ground . 

Un d c-rrro u nd  S t o r r g e . Storage  in  a cavern  or  chamber 
storage  site  with  a ceiling  more  than  60  cm  below  the 
natural  ground  level  and  storage  in  a site  tunnelled 
into  a hill. 


Und.oraro-.nd  Storage  Site  : A cavern  or  chamber  storage  site  v: 
sTceTlihs  not  less*  than  Go  cm  belo.v  the  natural  ground 
level , special!'-  adapted  for  the  storage  of  mnunition 
and  explosives. 


Venting  : The 

redact  io: 

n of  pressure 

complex 

due  to  t 

he  release  of 

another 

Q F fc  •‘••j  "r,  0 Y' 

through  an  a 

in  the  o 

vsrburde: 

.n  a chamber 
rases  into  a 
. t or  through 


or  storage 
corridor, 
or.  a Denture 


Vortex?.  A turbulent  movement  of  gas 
passing  constrictions,  corners, 


produced  by  blast  waves 
exits,  etc. 


Vise o u.s  log 
which 


s : A reduction  in  peak  pressure  at  the 
"results  from  the  roughness  of  tunnel  or 


hock 

tube 


front 
wails . 


V/?ll 


Ro  irnness  : The  average  height  of 

’oF~tiTe~we.il  of  a tube  or  tunnel. 


the  irregularities 


s>  cm  p> 
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CALIBRATION  OP  PRESSURE  TRANSDUCERS  AND  MEASURING  CHAIN 

The  purpose  of  the  calibration  was  to  check  the  pressure 
transducers  and  to  find  the  errors  introduced  by  the 
total  signal  processing  chain  and  digitizing  system. 

Pig.  B.1  shows  the  explosive  charge  (TNT)  driven  shock 
tube  used  in  the  calibration.  Six  pressure  transducers 
were  mounted  head  on  at  one  end  of  the  tube  on  a circle 
around  the  centreline  of  the  tube » The  transducers  were 
also  heat  protected  from  the  detonation  gases  by  a layer 
of  silicon  grease. 

The  measuring  chain  was  as  shown  in  Fig.  4.1  c (Mode  A), 
and  the  pressure-time  recordings  were  digitized  as  dis- 
cussed in  Appendix  C.  The  peak  pressure,  impulse  and 
positive  duration  as  well  as  mean  standard  deviations • 
were  found  employing  special  computer  programmes,  Appendix  C. 
The  calibration  data  supplied  by  the  manufacturer  were 
used  in  these  calculations.  Two  sets  of  six  transducers 
were  listed,  one  set  with  serial  number  and  higher  numbers. 

The  results  are  shown  in  Table  B.'i.  As  may  be  seen,  the 
mean  standard  deviation  for  peak  pressure  v/as  about  4# 
and  for  impulse  about  8$. 

During  the  course  of  the  experiments  to  be  discussed  in 
the  following  reports,  occasional  checks  of  this  type 
showed  that  the  transducer  calibrations  remained  remark- 
ably cons  tan.  t (-  5w) . There  were  a few  instances  where 
one  transducer  would  suddenly  fail,  but  this  could  easily 
be  detected  from  the  grossly  erroneous  pressure  data  and 
certain  characteristics  in  the  total  pressure- time  history 
(noise,  discharge  patters). 
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APPENDIX  C 


DATA  REDUCTION 

The  pressure-time  recordings  were  displayed  in  analog 
form  either  on  Polaroid,  plates  or  UV  sensitive  paper 
(Sec.  4.1),  which  were  not  suitable  for  further  analysis. 
The  analog  signals  were  therefore  converted  into  digital 
form  using  a D-mac  curve  tracer.  The  procedure  is  out- 
lined in  Pig.  C.1,  where  a relatively  small  number  of  dis 
•crete  points  on  the  curve  reproduces  the  essential  fea- 
tures of  the  pressure-time  history.  The  curve  tracer 
then  simply  converts  the  x,  y coordinates  for  each  of 
these  points  to  voltages,  which  are  digitised  using  an 
analog-to-digital  converter.  In  practice  up  to  approxi- 
mately 60  points  on  each  curve  digitised.  In  addition, 
a total  of  five  fixed  points  (1  - 5 in  Fig.  C.1)  were 
digitized  to  give  the  pressure  unit  (pts  1-2)  and  time 
unit  (pts  2-3).  Points  4 and  5 were  used  to  establish 
the  zero-pressure  line  (x-axis). 

The  digital  data  characterizing  each  recording,  were 
subsequently  processed  on  an  IBM  370/155  computer  and 
four  parameters  were  determined  employing  special  com- 
puter programmes: 


- Peak  pre ssure  (p) : The  maximum  value  of  the  pressure- 
time curve  obtained  by  taking  the  mean  values  of  the 
pressure  oscillating  system 

- Impulse  (I):  The  sum  of  trapezoids  between  the  x-axis 

and  pairs  of  successive  digitized  points  on  the  curve 
located  above  the  zero  over- pressure  line  (x-axis) 

- Positive  duration  (t , ) : The  time  difference  between  tb. 

arrival"' ’af~”£he ^shoek'Jront  (point  5 in  Pig.  C.1)  and 
the  crossing  of  the  smoothed  curve  with  the  zero  over- 
pressure line  (x-axis) 

- Arrival  time  (tp):  The  time  difference  between  the  re- 
cording's oi™tne  shock  fronts  in  the  detonation  chamber 
and  the  actual  measuring  point. 


C 2 


These  parameters  were  then  sorted  according  to  test  co 
figuration,  charge  weight  etc.  and  stored  for  further 
analysis . 


APPENDIX  D 


PITTING  TECHNIQUES  TO  ESTABLISH  SCALING  RELATIONSHIPS 

This  appendix  has  been  included  in  order  to  outline  a 
special  fitting  technique  used  in  Reports  III  and  IV,  to 
probe  various  empirical  scaling  relationships  and  obtain 
a consistent  set  of  scaling  parameters.  This  technique, 
originally  proposed  by  Powell  / 33/ , follows  the  normal 
least  squares  fitting  routine  up  to  certain  point. 


Let  x^,  y\  be  the  observed  data  to  be  fitted  to 
y = f(x,  a^ , b,  c,  ....),  where  f is  a function  described 
by  the  parameters  a,  b,  c,  ...  If  the  calculated  value 
of  y is  given  by  y.  - f (x, , a,  b,  c ...),  the  aim  of  the 
usual  least  squares  fit  is  to  minimise  L = (y.  - y ) . 


^_L 

a 


=)b  " 


= 0 


(D.1) 


In  the  least  squares  fitting  procedure,  these  equations 
are  solved  to  give  the  best  fit  values  of  a,  b,  c„  ... 
However,  the  function  f(x,  a,  b,  c,  .....)  may  be  so 
complex  that  Sc s. (D.1 ) cannot  be  solved  in  a straight- 
forward way. 

In  this  case  one  can  use  a method  which  in  essence 
approximates  the-  derivatives  by  differences  through 
an  iterative  procedure  /33/. 

In  practice  this  iterative  fitting  procedure  was  carried 
out  using  an  IBM  370/155  computer  and  the  programme  was 
written  to  be  able  to  fit  up  to  nine  parameters,  but 
actually  no  more  than,  four  were  used  in  the  present  -work, 
'The  programme  also  produced  standard  errors  on  the  fitted 
coefficients  as  well  as  their  correlations. 


APPENDIX  E 


E 1 


FILTER  TECHNIQUE  FOR  DETERMINATION  OF  MAXIMUM 
PRESSURE 


average 


This  appendix  describes  the  use  and  construction  of  two 
RC-filters  to  produce  the  maximum  value  of  a ■-  exponential 
decaying  signal  containing  severe  noise  or  ringing.  The 
technique  is  based  on  earlier  work  by  C.  Swensson  at  the 
Royal  Swedish  Fortification  Administration  / 32/,  but 
with  minor  changes  which  warrants  the  description  here. 

The  filter  technique  was  used  to  obtain  an  unambigous 
interpretation  of  the  pressure-time  recordings  in  the 
detonation  chamber  (Report  II).  These  recordings  were 
characterised  by  high  frequency  ringing  due  to  the  multip 
reflected  waves  from  the  chamber  walls  superimposed  on 
an  approximately  exponential  pressure  decay. 

The  high  frequency  ringing  can  be  .filtered  out,  but  then 
the  maximum  average  pressure  is  also  reduced.  As  will 
be  discussed,  it  is  possible  to  restore  the  true  average 
maximum  pressure  by  using  two  different  filterfrequences 
and  appropriate  correction  factors. 

Fig.  E.1  shov/s  the  principle  for  the  experimental  set-up. 
The  unfiltered  recorded  signal  is  fed  through  two-pole 
RC-filters  with  different  cut-off  frequencies 
f.  = 2 7/7 R.j  and  i?,  - 2 '/T/R^C^,  and  displayed  on  an 

UV  writer.  The  unfiltered  signal  is  also  displayed  for 
comparison. 


In  practice,  the  RC-filter  in  Fig.E 
filters  in  series  of  the  type  shown 
four  RC-filters  are  isolated  from  e 
operational  amplifier  of  AD  741  C1T 
follower  with  an  input  impedance  of 


. 1 consists  of  four- 
in  Fig.E.  2.  Tiie 
ach  other,  using  an 
the  type  as  a,  volts 
2 Mohrn. 


e 


The  following  components  were  chosen: 

Filter  no  \ 
j 1 

Vf  j 

x\i  ! 

(kJTl)  ! 

°3 

(A-f') 

-bud  freouencv= 

•J  u I 

(Hz ) 

1 

1C 

A 

1 

15,85 

2 

20 

1 

7,96 

With  four  filters  of  these  types  in  series  the  response 
function  has  a -12  dB/octave  high  frequency  roll-off. 

The  filters  were  applied  to  the  playback  signals  from 
the  Ampex  tape  recorder.  The  tape  speed  as  the  sig- 
nals were  recorded,  was  60  in  per  sec  and  the  play- 
back speed  was  1 7/8  in  per  sec.  This  means  a ..'actor 
of  '52  reduction  in  frequency  from  the  real  time  frequencie 
which  justifies  the  use  of  the  rather  low  filter  frequen- 
cies in  the  preceding  table. 

In  order  to  calculate  the  response  functions  at  a and  b • 
in  Fig.E.1  of  the  two  filtered  signals,  a Laplace  trans- 
formation is  used.  The  unfiltered  signal  is  assumed  to 
be  of  form 

p(t)  = PQe"  + n(t),  (3.1) 

where  n(t)  represents  the  superimposed  noise  as 
a function  of  time.  The  noise  is  assumed  to  be 
of  such  high  frequency  that  it  v/ill  be  removed 
by  the  tv/o  filters.  The  response  function  of  the 
four  series-coupled  circuits  in  Fig. 3. 2 is  then 


R(s)  = (, 


i=1 


¥i  • , 

— y"-)  s +T 


(3.2) 


O + T 


R.C. 

l .1 


where  R^Ch  is  the  time  constant  for  the  ith 
stage  (R1C1  = - R-C,.  - = RC). 

The  inverse  Laplace- transformation  produces  two  cases 
for  the  response  function  in  the  time  domain: 


L(t)  = po(iCTT) 

J d 


4 


2|  b e 


R.C 


•J  . 1 


(E.3) 


J 


1 or  2 corresponds  to  filter  1 or  2 in  Fig.E.1. 


E 3 


i 4 c 

P-j(t)  = P0(pf^Q™)  I e^P  (-t/'R.;C.) 


3 3 
t 


l 


(r-  i/r,c.)5  (« - i/a.c,)' 

J J JJ 

+ . exp  ( - cT/  t ) — j \ 

(eL-  1/R.C.)4J 


3 3 


1 


+3 


t2 


w 


1+R,CJ 
3 <> 


2 (C-  i/R.C, 
J J 


(E.4) 


If  A and  B are  the  maximum  values  of  the  filtered  signals 
in  Fig.E.I,  then  from  Eq.(E.3)  or  (E.4): 


A = p^t  =*C) 

E = p2(t  = <£)  (E.5) 

By  convention  we  shall  assume  A > B v/hich  is 
consistent  with  the  present  filter  frequencies. 
Clearly,  A < pQ  and  B < p and  it  is  convenient  to 
introduce  two  correction  factors  3c^  and  k9 
defined  as: 


from  which 


1 

1 


Wk„  - 1 

a = KA 

kg/K.j  ~ A/J5 


where 


a.  = 


i 


With  the  present  filter  frequencies, 
shows  K as  a function  A/B. 


(E.S) 


( E . 7 ) 


(E.S) 


3.3 


E 4 


As  an  illustration,  Fig.E.4  shows  the  .filtered 
signals  from  an  actual  experiment.  The  ratio 
"between  the  amplitude  of  the  two  filtered  signals 
is  found  to  be 

A/B  = 1,02 

Now , turning  to  Fig.E.3,  this  value  for  A/B 
produces  a correction  factor 

K = 1,025  ■ 


The  true  maximum  average  pressure  pQ  can  now  be  found 
from  Eq. (E, 7)  as : 

pQ  = IC  * A = 1,025  x 4,0  = 4,1  bar 

The  adequateness  of  this  whole  procedure  was  checked 
by  varying  the  filter  frequencies  f^  and  f2  over  a 
relatively  wide  range.  The  calculated  values  for 
appeared  to  be  reproducible  within  approximately 
5 - 10$.  An  increased  uncertainty  was  observed  for 
signals  which  deviated  from  the  ideal  exponential 
pressure  decay  in  Eq.(B.l), 

Essentially  the  same  procedure  can  be  followed  to 
estimate  the  positive  impulse 

I " ( Pn  s'  6/0  dt  = P J€  , 

-'G  U u 

but  this  wa.s  not  used  in  the  present  experiments. 
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Table  4.2.  System  specif ieat.ions 


Transdvicer  Range  Linearity 

(bar) 


Kistler 

6201 

0-5000 

+ 

0,4$ 

Kistler 

601 

H 

0-1000 

+ 

0,4$ 

Kistler 

603 

B 

0-  200 

+ 

0,5$ 

Charge  amplifiers 

Kistler  504  A,  565  and  568  Linearity  - 0,1  $ 
Kistler  5001  Linearity  - 0,05$ 

Tape  recorder  Ampex  PR  1300  A 
Voltage  linearity  - 0,75$ 

Harmonic  distortion  - 1,8$ 

Tape  speed  distortion  - 0,25'?° 

Galvanometer  amplifier 

Kistler  5211  Linearity  - 0,4$ 

U-V  writer 

CEO  type  124-5  Linearity  - 1 $ 


M-r 


time  xro:n  detonation  ox  charge  Q ^ 

TTR-J)  = loading  density  = 100  kg  Ti\:7//i;P 

/ 

/ 

/ 

/ 

/ 

/ 


Scale  Is  10  (Steel) 


Granite 


Lam£ ' s 

constant 

Elastic 

modulus 

V 

(N/m2) 

V 

(N/ra2) 

E 

(N/m2) 

•nr 

8,1  x 10,u 

12,1  x 1010 

2,1  x 1011 

EEBBB9HEEBBKI 

BEHBBH 

t 


a M.=  • ■■  , where  V = Poissons  ratio 

/ 1 + V 


E • V 

(1  + V ) (1  - 2v>  ) 


Fig.  2.3  Definition  of  symbols  and  numerical  values 
used  in  the  calculation  elastic  deforma- 
tion energy  for  a spherical  chamber. 
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V. 

1 

(cm^) 

A. 

D 

(cm2) 

Q 

4300 

^0ol 

TNT 

7250 

35.0 

PETR 

10900 

70.0 

AN/FO 

15200 

81.0 

RDX 

17300 

140.0 

ALUMIT 
DYNAMITE 
COMP.  B 

Fig. 3.1:  Model  configuration  used  in  the  closed  bomb 

tests  (Report  II),  with  different  chamber 

volumes  V. , vent  areas  A.  and  type  of  explosive  Q, 
1 J 


\Ac  =31 4 cm*’ 

V\\\W\ WWWWWNN  V 


\\v\\\\\\\\\ 


p 

V. 

vJ; 

1 

i 

(cn^) 

(mJ) 

7250 

330- 3100 

1C  900 

530-h7C0 

1 7200 

300-6500 

| 0 

A . 

3 

Co-,.2) 

2 

3 

(r2) 

1 

37 

20' 

2 

70 

20 

3 

InO 

20 

Fig. 3® 2:  Model  configurations  used  in  the 
storage  tests  (Report  III).  The  1 
volumes  V.  and  the  three  tunnel  cross 
A.  given  in  the  table  thus  io ruc- 
tions V.  A.  with  i - 1,  2,  3 and  ; 

1 3 2 

Choosing  K.  = 20  n as  a typical 
for  the  tunnel  crocs  section,  the 
factors  n have  the  values  shown  j 
column,  "he  corresponding  full  sc 

p 

volumes  owe  v . . 


Pig.  3.5  a Details  of  a typical  detonation  chamber. 


Copy  available  to  DDC  does  not 
permit  fully  legible  reproduction 


Pig.  3.5  c 


Details  of  tube  connections 


■it, 


.Copy  available  to  DDC  does  not 
permit  fully  legible  reproduction 


Detonat ion  c hamb er 


Silicon  oil 


Adapter  (steel) 


Kistler  pressure 
gauge 


Fig.  4.1  a Heat  protection  of  transducer  in  detonation 
chamber . 


v 71.1 

r-  ■■  " 

Jr 

! 

/ 

\ 

V 

U 

J 

A - — jU 

• >r  -• 

LJ 


Adapter 


2 mm  Silicon  rubber 


Fig.  4.1  b Section  of  tunnel  with  pressure  transducer. 


PRESSURE  TRANSDUCERS 


MEASURING  MODE  B 


Pig. kc 1C:  Measuring  chain  indicating  the  two  different 
modes  used  in  the'  tests. 


